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Abstract  o*  i - - - - - - - r 


In  a  previous  publication  "Hj;  we  described  the  power  spectral  densities 
(PSD's)  of  very  high  resolution  ( 1 0  m)  scalar  winds  in  the  stratosphere  based 
on  the  triangulated,  rocket-laid  smoke  trail  method.  The  PSD's  were  obtained 
there  by  means  of  a  fast  foutier  transform  version  of  the  Blackman-Tukey 
technique  pjftn  the  present  paper  we  shall  apply  the  same  methods  in  order 
to  obuin  the  PSD's  of  the  individual  velocity  components.  We  found  that  the 
slopes  and  amplitudes  of  the  spectra,  plotted  on  log-log  paper,  were  very 
similar  to  the  results  of  the  former  analysis  of  scalar  winds.  In  the  case  of  the 
/  latter,  u  obtained  the  average  slope  of  2.7  „£  0.2  and  log-amplitude  of 
''  3.24^0.2.  The  present  results  for  the  components  did  not  differ  in  a 

sutistically  significant  way  from  this.  On  the  other  hand,  we  found  a  statisti¬ 
cal  diffeience  between  the  amplitudes  of  zonal  and  meridional  component 


1.  Introduction 

In  a  previous  publication  [3]  the  power  spectral  densities 
(PSD’s)  were  presented  of  10  m  resolution  vertical  profiles  of 
horizontal  wind  speeds  in  the  stratosphere.  Figure  I  shows 
the  five  PSD’s  plotted  on  one  set  of  axes.  As  can  be  seen,  this 
figure  shows  the  universality  of  both  slope  and  amplitude  of 
these  PSD’s  very  strikingly.  The  velocity  fluctuations  are  here 
assumed  to  be  due  to  atmospheric  gravity  waves.  The  univer¬ 
sality  shown  in  Fig.  I  is  not  surprising  when  one  compares  it 
to  the  findings  made  in  the  upper  ocean.  In  that  medium  it  is 
well  known  that  the  gravity  waves  have  a  universal  spectral 
property  described  by  the  “G.M.”  or  “Garrett-Munk 
Model”  [4,  5].  The  latter  model  was  first  associated  with 
atmospheric  velocity  fluctuations  by  Dewan  [6].  It  was  how¬ 
ever,  Van  Zandt  [7]  who  first  proposed  that  the  atmosphere 
displayed  the  universal  aspects  of  this  model  and  he  has 
presented  a  significant  body  of  evidence  to  demonstrate  this. 

At  the  suggestion  of  R.  A.  Vincent  (private  communi¬ 
cation)  we  have,  in  the  present  paper,  applied  the  same  analy¬ 
sis  techniques  to  the  individual  velocity  components  as  were 
used  in  Ref.  [3]  in  connection  with  the  speeds  analyzed  there. 
The  main  question  we  address  here  is  whether  or  not  these 
velocity  component  PSD’s  differ  significantly  from  the  speed 
PSD’s.  As  will  be  shown,  the  answer  is  that  there  is  no 
statistically  significant  difference.  On  the  other  hand,  there  is 
evidence  for  a  difference  in  spectral  amplitudes  between  zonal 
and  meridional  components. 

2.  Experimciital  methods 

Our  data  were  obtained  from  rocket  laid  vertical  smoke  trails. 
These  trails  served  as  “fluid  markers”  and  they  were  in  the 
form  of  “dashes”,  the  height  ranges  of  which  are  given  in 
Table  I.  Time  lapse  photographs  taken  from  ground  based 
telescopic  cameras  were  used  to  triangulate  trail  positions  as 


10'* 


,0-»i - - - - - - - 1 

Cf»  KT*  icr’  ©■*  O'’ 

WAVE  NUMBER  (CY/m) 

SPEED  PSD'S 

Fig.  I.  Power  spectral  densities  (PSD’s)  of  speed  (scalar  wind)  profiles.  See 
Table  I  for  labels. 


Table  I .  Speed  profiles  (cubic  detrend)  PSD  parameters 
based  on  straight  One  fit  (from  k  =  10~^m~'  to  k  =  J/2 
Nyquist) 


Profile 

Height  range  (km) 

PSD  slope 

PSD  amplitude 
at  *  -  10'*c/m 
units:  (c/m) 

A  May  22 

(1)  29.3-33.4 

(2)  35.8-42.3 

25.9  ±  0.08 

1558 

B  May  20 

(1)  22.2-27.0 

(2)  29.5-33.5 

2.83  ±  0.05 

1970 

C  Apr  26 

(1)  18.4-21. 1 

(2)  23.5-27.8 

(3)  30-33.7 

2.74  ±  0.03 

1886 

D  May  2 

(1)  19.8-23.8 

(2) 28.0-31.3 

(3)  34.9-37.6 

2.93  ±  0.03 

1889 

E  Sep  12 

(1)  15.8-20.4 

(2)  23.1-26.9 

(3)  29.1-32.2 

(4)  34.1-36.6 

(5)  38.2-39.4 

2.46  ±  0.02 

1388 
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2.71  ±  0.187  1.738  ±  247 
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PROFILE:  "MAY  20* 


PROFILE:  *MAY  20* 


ICRinONAL  VELOQTY 


Fig.  2.  Example  of  velocity  component  profiles.  A.  Zonal;  B,  Meridional. 
This  is  Tor  the  “May  20”  profile. 


a  function  of  time  and  altitude.  These  in  turn  led  to  velocity 
profiles,  an  example  of  which  is  given  in  Fig.  2.  Further 
details  of  the  experimental  methods  were  published  in  Ref.  [3] 
and  thus  they  will  not  be  repeated  here.  In  short,  we  obtained 
what  is  considered  the  “state  of  the  art"  (i.e.,  10  m  resolution) 
velocity  profiles  derived  from  such  smoke  trails. 


3.  Method  of  data  analysis 

The  following  is  a  description  of  the  procedure  that  we  used 
to  spectrally  analyze  our  data.  It  is  identical  to  the  one  used 
for  Ref.  [3]. 

(1)  The  data  were  “prewhitened”  (Ref.  [8])  by  “first  dif¬ 
ferencing”  the  data  as  follows: 


(4)  A  “window”  or  “taper”  is  then  applied  (90%): 

Ar(Z)  =  ^(Z){l  -  +lj}  (zero  outside).  (5) 

(5)  The  “whitened”  (PSD.)  is  obtained  from 
PSD,(^)  =  2AZy'^^Ar(Z)cxpl2nj(kz/M)] 


I  ^  k  ^  MU. 

(6)  The  “post  darkened”  PSD  is,  finally, 
PSD.(A/3/AZ) 


PSDf-L.^  = 

\MAZ)  [2  - 


2  cos  (Inkfm)] 


(6) 

(7) 


K  =  w,  -  o,., 


4.  Results 


where  v,  are  the  velocity  component  values,  /  is  the  data  point 
number,  and  y,  represent  the  prewhitened  data. 

(2)  The  mean,  trend,  and  rate  of  change  of  trend  were  then 
removed  by  subtracting  out  a  fitted  cubic  polymomial. 

(3)  The  autocorrelation  was  obtained  from  the  n  points  by 
first  adding  zeros  until  there  were  M  points  such  that 
M  =  2y  >  2N.  Ther  y  is  an  integer  and  this  is  done  so  that 
we  can  apply  a  Fast  Fourier  Transform  (FFT).  Next  we  used 
the  following  FFT  technique  to  obtain  the  autocorrelation. 
First  we  computed 

=  Yv,txp[2nAki)IM]  (2) 

where  AZ  is  the  altitude  separation  between  points  (10  m 
most  of  the  time).  The  autocorrelation,  which  is  A(Z),  is 
defined  by 


Figure  3  summarizes  the  PSD’s  of  the  meridional  and  zonal 
profiles.  As  can  be  seen.  Figs  I  and  3  (A,  B)  resemble  each 
other.  Table  II  (A,  B)  displays  the  individual  slopes  and 
amplitudes  of  the  least  square  fitted  straight  lines  extending 
from  wavenumber  k  =  10"*m''  to  k  =  1/2  Nyquist 
wavenumber,  as  well  as  to  =  (1/100  m)  (compare  Table 
I  for  the  results  for  the  speeds).  Tables  I  and  II  also  give  the 
averages  and  standard  deviations  of  the  slopes  and  ampli¬ 
tudes. 

A  “r-test”  [10]  was  applied  to  these  averages  to  see  if  there 
was  a  significant  difference  between  the  (a)  meridional  and 
zonal  component  slopes  and  amplitudes  (at  k  -  I0~'m~') 
and  (b)  the  velocity  component  slopes  and  amplitudes  vs. 
those  of  the  speeds.  Table  III  summarizes  the  results.  To 
apply  the  r-test  we  used,  as  estimate  of  the  error  of  the 
means,  a. 


(tt  jjj .  r ("i  -  i)<t?  +  ("j  -  i)<H] 

L^J  y-o 

L  rii  Hj  —  2  J 

L"!  "ij 

A(Z) 

is  obtained  from  (see  Ref.  [9]) 


A{x) 


A/AZ 


exp  [-2nj(kzlM)] 


'}■ 

(4) 


where  a^  and  Oj  were  calculated  from  the  samples  from 


n  —  1 


(8) 


(9) 
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Fig.  3.  Power  spectral  densities  (PSD's)  of  Zonal  (A),  and  meridional  (B) 
velocit>  components.  Table  1  shows  labels,  as  in  Fig.  I. 


Table  II.A.  Velocity  profiles  (cubic  detrend)  PSD  parameters  based  on  straight  line  fit  (from  k  =  10'^ m~'  to  k  =  112 
Nyquist) 


Profile 

Meridional 

Zonal 

PSD  slope 

PSD  amplitude 
at*  -  lO-’m*' 
units:  (m’/s')  (c/m) 

PSD  slope 

PSD  amplitude 
at  *  “  I0‘’m"’ 
units:  (m'/s*)  (c/m) 

May  22 

2.S9 

1069 

2.78 

2402 

May  20 

2.57 

995.4 

2.99 

2243 

Apr  26 

2.73 

1951 

2.65 

1764 

May  2 

3.06 

1618 

2.93 

1779 

Sep  12 

2.43 

1642 

2.53 

1645 

Averages 

2.68  ±  0.24 

1455  ±  409 

2.78  ±  0.19 

1967  ±  334 

Combined  averages: 

2.73  ±  0.21 

1711  ±  443 

(both  directions 

Table  II.B.  Velocity  profiles  (cubic  detrettd)  PSD  parameters  based  on  straight  line  fit  (from  k  *  10~^m~' 

'  to  k  ^  1/lOOmJ 

Profile 

Meridional 

Zonal 

PSD  slope 

PSD  amplitude 

PSD  slope 

PSD  amplitude 

at  t  -  IO-»m-' 

at  *  “  I0"’m'‘ 

units:  (m’/s*)  (c/m) 

uniu:  (m’/s^)  (c/m) 

May  22 

2.65 

1141 

2.94 

2855 

May  20 

3.03 

2028 

3.22 

3110 

Apr  26 

2.79 

2072 

3.00 

2791 

May 

3.05 

1623 

3.52 

4218 

Sept  12 

2.63 

2183 

2.88 

2503 

Avenges 

2.83  ±  0.202 

1809  ±  430 

311  ±  0.262 

3095  ±  664 

Combined  averages 

2.97  +  0.266 

7457  +  859 
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Table  III.  t-Test  between  means  of  PSD  parameters  ("significance”  means  ^5%  level  of  significance) 


Type  of  proiile 

Slope 

Amplitude 

Value  of  r 

Sig.  level 

Value  of  T 

Sig.  level 

Meridional  vs.  zonal 

0.73 

50% 

2.17 

10% 

(-U.  -  NYQ/2) 

df-  8 

Not  sig. 

df  -  8 

Not  sig. 

Meridional  vs.  Zonal 

1.89 

10% 

3.64 

1% 

(2»  -  100m) 

df  •  8 

Not  sig. 

df  -  8 

Signif. 

Vector  components  vs.  Speed 

0.180 

90% 

0.125 

90% 

(-!«.  -  NYQ/2) 

df  <=  13 

Not  sig. 

df  -  13 

Not  sig. 

Thus  for  t  we  have: 

_ -  2(P,  -  P,) _ 

y/OM  +  (iM  v/x  yf  -  n,(p,9  +  ^y}-  nM' 

(10) 

As  can  be  shown  from  and  average  variances,  there  are 
no  significant  difierences  between  component  and  speed  aver¬ 
age  a^.  Note  that,  consistent  with  this  finding,  the  PSD’s 
(Table  III)  have  no  statistically  different  parameters  from  one 
another  with  one  exception.  The  zonal  and  meridional  spec¬ 
tral  amplitudes  do  differ  statistically,  when  we  take  the  cut  off 
wavelength  as  100  m.  This  is  an  imporunt  result  in  the  sense 
that  universality  does  not  apply  to  the  component  spectra  but 
rather  to  the  total  spectra.  Theoretically,  this  was  to  be 
expected. 

6.  Conclusion 

High  resolution  vector  wind  profiles  in  the  stratosphere  have 
been  spectrally  analyzed  and  compared  to  speed  profile  spec¬ 
tra  previously  analyzed  in  Ref.  [3].  The  meridonal,  zonal  and 
speed  PSD’s  were  compared  in  regard  to  the  slope  and  the 


amplitude  (at  k  »  10~’m~')  and  no  significant  differences 
were  found  between  speeds  and  velocity  components.  There 
was  however  statistical  evidence  for  an  amplitude  difference 
between  zonal  and  meridional  component  PSD’s.  Our  main 
conclusion  is  that  there  is  no  statistical  difference  between  the 
power  spectral  densities  of  velocity  components  taken  as  a 
whole  and  the  PSD’s  of  speeds. 
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